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I. INTRODUCTION 

Double-logarithmic (DL) contributions are of a special interest among radiative corrections. They are interesting 
in two aspects: first, in every fixed order of the perturbation theories they are the largest terms among the radiative 
II . corrections depending on the total energy and second, they are easiest kind of the corrections to sum up. DL corrections 
Q_i' were discovered by V.V. Sudakov in Ref. [l| in the QED context. He showed that DL terms appear from integrations 
over soft, infrared (IR) -divergent momenta of virtual photons. All-order resummation of such contributions led to 
their exponentiations. 

Next important step was done in Refs. [2[ where calculation and summation of DL contributions was considered 
in a systematic way. They found a complementary source of DL terms: soft virtual fcrmions. This situation appears 
in the Regge kinematics. The all-order resummations of DL contributions in the Regge kinematic are quite involved 
and yield more complicated expressions than the Sudakov exponentials. Nonetheless important was the proof of the 
factorization of bremsstrahlung photons with small k± in the high-energy hadronic reactions found in Ref. [3| and 
often addressed as the Gribov's bremsstrahlung theorem. This statement, suggested originally in the framework of 
the phenomenological QED of hadrons was extended to QCD in Refs. [3]. 

Calculation in the double-logarithmic approximation (DLA) amplitudes of the fermion-antifermion annihilation in 
the Regge forward and backward kinematics involves accounting for DL contributions from soft quarks and soft gluons. 
These reactions in QED and QCD have many common features. The e + e~ -annihilation was studied in Refs. 0. 
The quark-aniquark annihilation DLA was investigated in Ref. The method of calculation here was based on 
, ' factorization of virtual quarks and gluons with minimal k±. Generally speaking, the results obtained in Ref. [fj( could 
^ . be obtained with the method of Ref. 0, however the technique of calculations suggested in Ref. @ was much more 
elegant and efficient. Although Ref. [5[ is about quark scattering only, it contains almost all technical ingredients 
necessary to compose Infrared Evolution Equations for any of elastic scattering amplitudes. Nevertheless it could not 
directly be applied to inelastic processes involving emission of soft particles. Such a generalization was obtained in 
Refs. @, Q. The basic idea of the above-mentioned method was suggested by L.N. Lipatov: to investigate evolution 
with respect to the infrared cut-off. The present, sounding naturally term "Infrared Evolution Equations" (IREE) for 
this method was suggested by M. Krawczyk in Ref. Q where amplitudes for the backward Compton scattering were 
calculated in DLA. 

The aim of the present brief review is to show how to compose and solve IREE for scattering amplitudes in different 
field theories and kinematic regions. The paper is organized as follows: in Sect. II we consider composing IREE in the 
technically simplest hard kinematics. In Sect. Ill we consider composing IREE in the forward kinematics and apply 
it to studying the structure function gi of the polarized Deep-Inelastic scattering (DIS) at small x. The point is that 
the commonly used theoretical instrument to study <?i is DGLAP [11]. It collects logarithms of Q 2 to all orders in 
a s but does not include the total resummation of logarithms of l/x, though it is important at small x. Accounting 
for such a resummaton leads to the steep rise of g\ at the small- a; region. As is shown in Sect. IV, DGLAP lacks the 
resummaion but mimics it inexplicitly, through the special choice of fits for the initial parton densities. Invoking such 
peculiar fits together with DGLAP to describe g\ at x <C 1 led to various misconceptions in the literature. They are 
enlisted and corrected in Sect. V. The total resummaion of the leading logarithms is essential in the region of small x. 
In the opposite region of large x, DGLAP is quite efficient. It is attractive to combine the resummation with DGLAP. 
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The manual for doing it is given in Sect. VI. Finally, Sect. VII is for concluding remarks. 



II. IREE FOR SCATTERING AMPLITUDES IN THE HARD KINEMATICS 



From the technical point of view, the hard kinematics, where all invariants are of the same order, is the easiest for 
analysis. For the simplest, 2 — > 2 -processes, the hard kinematics means that the Mandelstamm variables s,t,u obey 

s 1 u . (1) 

In other words, the cmf scattering angles 9 ~ 1 in the hard kinematics. This kinematics is the easiest because the 
ladder Feynman graphs do not yield DL contributions here and usually the total resummation of DL contributions 
leads to multiplying the Born amplitude by exponentials decreasing with the total energy. Let us begin with composing 
and solving an IREE for the well-known object: electromagnetic vertex Y „ of an elementary fermion (lepton or quark). 
As is known, 

r» = u(p 2 )[j,f(q 2 )-^ 9 (q 2 )]u( Pl ) (2) 

where pi ; 2 are the initial and final momenta of the fermion, m stands for the fermion mass and the transfer momentum 
q = p 2 — Pi- Scalar functions / and g in Eq. ^ are called form factors. Historically, DL contributions were discovered 
by V. Sudakov when he studied the QED radiative corrections to the form factor / at \q 2 \ ^> \p\ 2 |. Following him, 
let us consider vertex at 

\q 2 \^p\=p\=m 2 (3) 

i.e. we assume the fermion to be on-shell and account for DL electromagnetic contributions. We will drop m for the 
sake of simplicity. 



A. IREE for the form factor f(q 2 ) in QED 

Step 1 is to introduce the infrared cut-off \i in the transverse (with respect to the plane formed by momenta Pi^) 
momentum space for all virtual momenta ki : 

h ± > \i (4) 

where i = 1,2, ... 

Step 2 is to look for the softest virtual particle among soft external and virtual particles. The only option we have 
is the softest virtual photon. Let denote its transverse momenta = k±. By definition, 

fc_L = min ki ± . (5) 

Step 3: According to the Gribov theorem, the propagator of the softest photon can be factorized (i.e. it is attached 
to the external lines in all possible ways) whereas k± acts as a new cut-off for other integrations. Adding the Born 
contribution f Born = 1 we arrive at the IREE for / in the diagrammatic form. It is depicted in Fig. 1. IREE in the 
analytic form are written in the gauge-invariant way, but their diagrammatical writing depends on the gauge. In the 
present paper we use the Feynman gauge. 

Applying to it the standard Feynman rules, we write it in the analytic form: 

fin 2 - f B °™ e2 f dadpdkj e(fcj-^ 2 ) f(q 2 ,kp 



8tt 2 J (saj3 — k\ + ie)(—sa + saj3 — k\ + ie)(s[3 + sa(3 — k 2 _ + le) 

where we have used the Sudakov parametrization k = ap2 + Ppi + k± and denoted s = —q 2 w 2pip 2 . As f(q 2 , kj_) 
does not depend on a and (3, the DL integration over them can be done with the standard way, so we are left with a 
simple integral equation to solve: 



/(iV) = /""*' -n ? ln(»Al)/(9 2 , *i) . (7) 
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FIG. 1: The IREE for the Sudakov form factor. The letters in the blobs stand for IR cut-off. 

Differentiation of Eq. ([7|) over p, 2 (more exactly, applying —p 2 d/dp 2 ) reduces it to a differential equation 

df/d(ln( S /p 2 )) = ~(e 2 /8n 2 ) ln( S //x 2 )/ (8) 

with the obvious solution 

/ = / So ™exp[-(a/4vr) \n 2 (q 2 /m 2 )} (9) 
where we have replaced /ibym and used a — e 2 /4n. Eq. ^ is the famous Sudakov exponential obtained in Ref. 

B. IREE for the form factor g(q 2 ) in QED 

Repeating the same steps (see Ref. Q for detail) leads to a similar IREE for the form factor g: 

g(q 2 ,m 2 ,p 2 ) = g Born (s,m 2 ) - * f §L ln( S /fc 2 ). 9 (g 2 , m 2 , fc 2 ) (10) 

where g Born (s,m 2 ) — —(m 2 /s)(a/ir)ln(s/m 2 ). Solving this equation and putting p = m in the answer leads to the 
following relation between form factors / and g: 

9(s) = ~2 d / p , (11) 
with p — s/m 2 . Combining Eqs. (|9I11[) allows to write a simple expression for the DL asymptotics of the vertex T^: 

T M = u(p 2 ) [ 7/1 + ^±^-] u ( Pl ) exph(a/47r) In 2 p] . (12) 
m op 

C. e + e~ -annihilation into a quark-antiquark pair 

Let us consider the e + e~ -annihilation into a quark q(p±) and ?(pa) & t high energy when 2p\p 2 3> p 2 2 - We consider 
the channel where the e + e~ -pair annihilates into one heavy photon which decays into the q(pi) q(j>2) -pair: 

e+e~ -> 7* -> g(pi) q(p 2 ) ■ (13) 

We call this process elastic. In this case the most sizable radiative corrections arise from the graphs where the quark 
and antiquark exchange with gluons and these graphs look absolutely similar to the graphs for the electromagnetic 
vertex considered in the previous subsection. As a result, accounting for the QCD radiative corrections in DLA 
to the elastic form factors f q , g q of quarks can be obtained directly from Eqs. (|9llip by replacement 

a -> a s C F , (14) 

with C F = (N 2 - 1)/2N = 4/3. 
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D. e e -annihilation into a quark-antiquark pair and gluons 



In addition to the clastic annihilation (|13|) . the final state can include gluons: 

e + e~ -> 7* -> g(p 2 ) + sOi), ••fl'(fcn) ■ (15) 



We call this process the inelastic annihilation. The QED radiative corrections to the inelastic annihilation (]T5f) in 
DLA are absolutely the same as the corrections to the elastic annihilation. On the contrary, the QCD corrections 
account for gluon exchangesbetween all final particles. This makes composing the IREE for the inelastic annihilation 
be more involved (see Ref. Q). The difference to the considered elastic case appears at Step 2: look for the softest 
virtual particle among soft external and virtual particles. Indeed, now the softest particle can be both a virtual gluon 
and an emitted gluon. For the sake of simplicity let us discuss the 3-particle final state, i.e. the process 



q(pi) q(p2)+9(h) . (16) 
i proces: 

In DLA, it is parameterized by new form factors F^ 1 ' and 



The main ingredient of the scattering amplitude of this process is the new electromagnetic vertex F^ 1 ' of the quark 



IV = £i(fci)£(p 2 ) [ 7 ^ (1) (9, fci) - ^G«(<z, fci)]«(p!) (17) 

where (1) corresponds to the number of emitted gluons, q = pi + P2 and I is the polarization vector of the emitted 
gluon. The bremsstrahlung factor B\ in Eq. (JT7J) at high energies is expressed through k\ j_ : 

\p 2 ki pikiJ ku_ 

We call fWiG'"' inelastic form factors. Let us start composing the IREE for F^\ Step 1 is the same like in 
the previous case. Step 2 opens more options. Let us first choose the softest gluon among virtual gluons and denote 
its transverse momentum k± The integration over k± runs from /j to s. As /i < k\ j_ < s, we have two regions to 
consider: Region D\ were 

H < k 1± < k± < (19) 

and Region D 2 were 

/i < kj_ < ki± < y/s (20) 

Obviously, the softest particle in Region D\ is the emitted gluon, so it can be factorized as depicted in graphs (b,b') 
of Fig. 2. 

On the contrary, the virtual gluon is the softest in Region D 2 were its propagator is factorized as shown in graphs 
(c,d,d') of Fig. 2. Adding the Born contribution (graphs (a, a') in Fig. 2) completes the IREE for depicted in 
Fig. 2. Graphs (a-b') do not depend on n and vanish when differentiated with respect to \i. Blobs in graphs (c-d') do 
not depend on the longitudinal Sudakov variables, so integrations over a, (3 can be done like in the first loop. After 
that the differential IREE for F^ is 



dfi 2 2tt 

Solving Eq. (f2Tj) and using that (2pifci)(2p2&i) = skf ± leads to the expression 

,2 

,11" J 2 \ [I 



suggested in Ref. 0] and proved in Ref. [i| for any n. The IREE for the form factor G (n) was obtained and solved in 
Ref. M. It was shown that 



G*(«) = -2BF^/dp 



(23) 
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FIG. 2: The TREE for the inelastic quark form factor. 



E. Exponentiation of Sudakov electroweak double-logarithmic contributions 

The IREE -method was applied in Ref. [10] to prove exponentiation of DL correction to the electroweak (EW) 
reactions in the hard kinematics. There is an essential technical difference between the theories with the exact gauge 
symmetry (QED and QCD) and the EW interactions theory with the broken SU(2) U(l) gauge symmetry: only 
DL contributions from virtual photons yield IR singularities needed to be regulated with the cut-off /x whereas DL 
contributions involving W and Z -bosons are IR stable because the boson masses M\y and Mz act as IR regulators. 
In Ref. [l(| the difference between M w and M z was neglected and the parameter 



M > M- 



w 



M, 



(24) 



was introduced, in addition to ^, as the second IR cut-off. It allowed to drop masses M\y,z- The IREE with two IR 
cut-offs was composed quite similarly to Eq. ([6]), with factorizing one by one the softest virtual photon, Z-boson and 
VF-boson. As a result the EW Sudakov form factor Few is 



F, 



ew = exp 



a{Q\ + Q 
8tt 



9 2 C F 
16tt 2 



SU(2) 



g' 2 (F 1 2 +y 2 2 ) a{Q\ + Ql) 



16tt 2 



Stt 



In 2 (s/M 2 



(25) 
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where Q12 are the electric charges of the initial and final fermion (with W -exchanges accounted, they may be 

different), Yi, 2 are their hyper-charges and Cp U{2) = (N 2 - 1)/2N, with N = 2. We have used in Eq. J35J the 
standard notations g and g' for the SU(2) and U(l) -EW couplings. The structure of the exponent in Eq. (f2"5|) is quite 
clear: the first, /i -dependent term comes from the factorization of soft photons like the exponent in Eq. ([9]) while 
other terms correspond to the W and Z -factorization; the factor in the squared brackets is the sum of the SU(2) 
and U(l) Casimirs, with the photon Casimir being subtracted to avoid the double counting. In the limit /i = M the 
group factor in the exponent is just the Casimir of SU(2) ® U(l). 



III. APPLICATION OF IREE TO THE POLARIZED DEEP-INELASTIC SCATTERING 

Cross-sections of the polarized DIS are described by the structure functions 31,2- They appear from the standard 
parametrization of the spin-dependent part W pv of the hadronic tensor: 



W,t V = te^uxpqx— S p gi_(x,Q 2 ) + [S p - p p — 102(2:, Q 2 ) (26) 
pql V pq J J 



where p, m and S are the momentum, mass and spin of the incoming hadron; q is the virtual photon momentum; 
Q 2 = -q 2 ; x = Q 2 /2pq. Obviously, Q 2 ^ and < x s$ 1. 

Unfortunately, g\^ cannot be calculated in a straightforward model-independent way because it would involve QCD 
at long distances. To avoid this problem, W pl , is regarded as a convolution of <& q , g - probabilities to find a polarized 

quark or gluon and the partonic tensors W^l' 9 ^ parameterized identically to Eq. (f2"6")l . In this approach W p t' 9 ^ involve 
only QCD at short distances, i.e. the Perturbative QCD while long-distance effects are accumulated in & q , g - As <fr q ^ g 
are unknown, they are mimicked by the initial quark and gluon densities 8q, 5g. They are fixed aposteriori from 
phcnomenological considerations. So, the standard description of DIS is: 

* W$ ®5q + W$ ® Sg . (27) 

The standard theoretical instrument to calculate gi is DGLAP[ll[ complemented with standard fits[l2| for 5q, 5g. 
We call it Standard Approach (SA). In this approach 



9i {x, 1 



C q (x/z) ® Aq(z, Q 2 ) + C g (x/z) ® Ag(z, 1 



(28) 



where Cq, g are coefficient functions and Aq(z, Q 2 ), Ag(z, Q 2 ) are called the evolved (with respect to Q 2 )quark and 

us 

P gq Aq + P gg Ag] (29) 



gluon distributions. They are found as solutions to DGLAP evolution equations 

dAq a s (Q 2 ) , dAg a s (Q 2 ) 

[P qq Aq + P qg Ag\ , 



d\nQ 2 



2tt 



d\nQ 2 



2tt 



where P a b are the splitting functions. The Mellin transforms 7 a t of P a & are called the DGLAP anomalous dimensions. 
They are known in the leading order (LO) where they are ~ a s and in the next-to- leading order (NLO), i.e. ~ a 2 . 
Similarly, C q . g are known in LO and NLO. Details on this topic can be found in the literature (e.g. see a review 
[HI]). Structure function gi has the flavor singlet and non-singlet components, gf and g^ s . Expressions for g^ s are 
simpler, so we will use mostly them in the present paper when possible. It is convenient to write g\ in the form of the 
Mellin integral. In particular, 
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NS DGLAP 1 



(x,Q 2 ) = (e 2 q /2) j" ^(~yC NS (cj)Sq(cj)exp 



Q dk 2 

-rrlNS^,a s {ki)) 

2 K 'I 



(30) 



where /i 2 is the starting point of the Q 2 -evolution; Cms and ^ns are the non-singlet coefficient function and anomalous 
dimension. In LO 



7ivs(w,Q 2 ) 



2tt 



Cff s (u) = l + 



a s {Q 2 )C F 
2vr 



1 

2uj 



1 



2uj + 1 



Lw(l + u>) 2 
(&(w)+^(w)-5 2 (w)) 



(31) 



w(l + u) 




+ 




FIG. 3: The IREE for the non-singlet component of the spin structure function g 1 . 

with S r (uj) = ^2j—i l/j r • The initial quark and gluon densities in Eq. (|30[) are defined through fitting experimental 
data. For example, the fit for 5q taken from the first paper in Ref. [l2j is 

5q(x) = Nx~ a \{1 - xf (1 + jx s )] , (32) 



with N being the normalization, a — 0.576, (3 — 2.67, 7 = 34.36 and 6 
DGLAP equations were suggested for describing DIS in the region 



0.75. 



x< 1, 



f » V 2 



(33) 



(fj, stands for a mass scale, fj, ^> Aqcd) and there is absolutely no theoretical grounds to apply them in the small- a; 
region, however being complemented with the standard fits they are commonly used at small x. It is known that SA 
provide a good agreement with available experimental data but the price is invoking a good deal of phcnomcnological 
parameters. The point is that DGLAP, summing up leading ln fe Q 2 to all orders in a s , cannot do the same with 
leading ln fc (l/x). The later is not important in the region (f33|) where ln fe (l/a;) -C 1 but becomes a serious drawback 
of the method at small x. The total resummation of DL contributions to <?i in the region 



x -C 1, 



(34) 



was done in Refs. [14j . The weakest point in those papers was keeping a s as a parameter, i.e. fixed at an unknown 
scale. Accounting for the most important part of single-logarithmic contributions, including the running coupling 
effects were done in Refs. [ll|. In these papers /1 2 was treated as the starting point of the Q 2 -evolution and as the IR 
cut-off at the same time. The structure function gi was calculated with composing and solving IREE in the following 
way. 

It is convenient to compose IREE not for g\ but for forward (with |t| < /i 2 ) Compton amplitude M related to gi 
as follows: 



91 



1 

7T 



(35) 



It is also convenient to use for amplitude M the asymptotic form of the Sommerfeld- Watson transform: 

27TI \fl A J 



M 



(36) 



where S}-^{uj) — [e~ llrw — l]/2 w —nruj/2 is the signature factor. The transform of Eq. (f3"6"]) and is often addressed 
as the Mellin transform but one should remember that it coincides with the Mellin transform only partly. IREE for 
Mcllin amplitudes F(oj, Q 2 ) look quite simple. 

For example, the IREE for the non-singlet Mellin amplitude F NS related to g^ s by Eqs. (|35l36p is depicted in 
Fig. 3. In the Mellin space it takes the simple form: 



w + d/dy]F Nb = (1 + lu/2)H ns F 



NS 



(37) 
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where y = ln(Q 2 //i 2 ). Eq. (|37|) involves a new object (the lowest blob in the last term in Fig. 3): the non-singlet 
anomalous dimension Hj^g accounting for the total resummaton of leading logarithms of 1/x. Like in DGLAP, the 
anomalous dimension does not depend on Q 2 but, in contrast to DGLAP, H^s can be found with the same method. 
The IREE for it is algebraic: 



coH NS = A(uj)C f /8tt 2 + (1 + lu/2)H 2 ns + D(lu)/%-k 2 



(38) 



The system of Eqs. (|37I38[) can be easily solved but before doing it let us comment on them. The left-hand sides 
of Eqs. (|37l38p are obtained with applying the operator — p 2 d/dp, 2 to Eq. (j36)) . The Born contribution in Fig. 3 
does not depend on p and therefore vanishes. The last term in Fig. 3 (the rhs of Eq. (|37[) ) is the result of a new, t 
-channel factorization which does not exist in the hard kinematics defined in Eq. {!]). In order to compose the IREE 
for the Compton amplitude M, in accordance with the prescription in the previous section we should first introduce 
the cut-off p. Then Step 2 is to tag the softest particles. In the case under discussion we do not have soft external 
particles. Had the softest particle been a gluon, it could be factorized in the same way like in Sect. II. However, the 
only option now is to attach the softest propagator to the external quark lines and get ln(i//i 2 ) = from integration 
over (3 (cf Eq. (O). So, the softest gluon does not yield DL contributions. The other option is to find a softest quark. 
The softest t -channel quark pair factorizes amplitude M into two amplitudes (the last term in Fig. 3) and yield DL 
contributions. The IREE for H^g is different: 

(i) H N s does not depend on Q 2 , so there is not a derivative in the lhs of Eq. (|3T|) . 

(ii) The Born term depends on p and contributes to the IREE (term A in Eq. (j3"7) )). 

(iii) As all external particles now are quarks, the softest virtual particle can be both a quark and gluon. The case 
when it is the t -channel quark pair, corresponds to the quadratic term in the rhs of Eq. (|37|) . The case of the softest 
gluon yields the term D, with 



D(u) 



2Cp r 

b 2 N J 



p + r) 



(p + rf) 2 + 7r 2 rj 



(39) 



where b = (33 - 2n/)/127r and rj = \n{p 2 / K 2 QC D ) . 

The term A in Eq. (|37[) stands instead of a s . The point is that the standard parametrization a s = a s (Q 2 ) cannot 
be used at x <C 1 and should be changed (see Ref. [16| for detail). It leads to the replacement a s by 



A( U ) = - 



dpe 



L?7 z +7r z J (p + r/) 2 +7r 2 
Having solved Eqs. (|37l38p . we arrive at the following expression for s in the region (|34 



/too J 
— (1/zrCjvsM^Hexp (H ns {lo) V ) 



where the coefficient function Cjvs(w) is expressed through Hns{u)'- 



u - H ns (uj) 



and Hns(^) is the solution of algebraic equation 



(40) 



(41) 



(42) 



H NS = (1/2) |w - ^Ju 2 - B{lo) 

where 

B(uj) = (47rC F (l + uj/2)A{uo) + L>(w))/(2tt 2 ) . 
It is shown in Ref. [17| that the expression for gi in the region 



(43) 
(44) 



x < 1, 



< 



/'" 



can be obtained from the expressions obtained in Refs. for g\ in region (f3~4"|) by the shift 



> 2 +/io 



(45) 
(46) 



where //q = 1 GeV for the non-singlet g\ and //q = 5.5 GeV for the singlet. 
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IV. COMPARISON OF EXPRESSIONS ([30]) AND (gT]) FOR ff f s 

Eqs. ([30]) and (j4"Tj) read that the non-singlet gi is obtained from Sq with evolving it with respect to x (using the 
coefficient function) and with respect to Q 2 (using the anomalous dimension). Numerical comparison of Eqs. (|30|) 
and l|4ip can be done when Sq is specified. 



A. Comparison of small-a; asymptotics, neglecting the impact of 5q 

In the first place let us compare the small- a; asymptotics of for gf s DGLAP and g^ s , assuming that Sq does not 
affect them. In other words, we compare the differencee in the rc-evolution at x — » 0. Applying the saddle-point 
method to Eqs. (|30[) and (|4ip leads to the following expressions: 



and 



NS DGLAP 

9i 



9? S 



exp 



^ln(l/x)lnln(QVA| CD ) 



(l/a;) Ajvs (Q 2 /M 2 ) Ajvs/2 



(47) 



(48) 



where A^s = 0.42 is the non-singlet intercept 1 . Expression (|4"T)) is the well-known DGLAP asymptotics. Obviously, 
the asymptotics is much steeper than the DGLAP asymptotics ([50)1 . 



B. Numerical comparison between Eqs. (|30p and (I41|) . neglecting the impact of <5q 

A comparison between Eqs. (|30p and (|41[) strongly depends on the choice of <5q but also depends on the difference 
between the coefficient functions and anomalous dimensions. To clarify the latter we choose the simplest form of Sq: 

Sq(uj) = N q . (49) 

It corresponds to the evolution from the bare quark where Sq(x) = N q S(l — n 2 /s). Numerical results for R = 
[9i S — 9i S DGLAP ]/g^ s dglap witn fig chosen by Eq. (J49]) manifest (see Ref. for detail) that R increases when 
x is decreases. In particular, R > 0.3 at x < 0.05. It means that the total resummation of leading ln fe (l/x) cannot 
be neglected at x < 0.05 and DGLAP cannot be used beyond x w 0.05. On the other hand, it is well-known that 
Standard Approach based on DGLAP works well at x <§C 0.05. To solve this puzzle, we have to consider the standard 
fit for Sq in more detail. 



C. Analysis of the standard fits for Sq 



There are known different fits for Sq. We consider the fit of Eq. (|32|) . Obviously, in the u> -space Eq. (|32|) is a sum 
of pole contributions: 

5q(w) = N V [(w - a)" 1 + £ m k (u + Afe)" 1 ] , (50) 

with Afc > 0, so that the first term in Eq. ([ST))) corresponds to the singular term x~ a of Eq. (|3"2"|) and therefore the 
small- a; asymptotics of /dglap is given by the leading singularity u = a = 0.57 of the integrand in Eq. (|50[) so that 
the asymptotics of gf s DGLAP (x, Q 2 ) is not given by the classic exponential of Eq. (|47| but actually is the Regge-like: 



NS DGLAP 

9i 



C( a )(l/x)«(ln(Q 2 /A 2 )/ln( M 2 /A 2 )) 7(a)/b , (51) 



with b = (33 — 2n/)/127r. Comparison of Eq. (|4^|) and Eq. (pjTj) demonstrates that both DGLAP and our approach 
lead to the Regge behavior of gi, though the DGLAP prediction is more singular than ours. Then, they predict 



1 The singlet intercept is much greater: Ag = 0.86. 
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different Q 2 -behavior. However, it is important that our intercept Ajvs is obtained by the total resummation of the 
leading logarithmic contributions and without assuming singular fits for Sq whereas the SA intercept a in Eq. (|47[) 
is generated by the phcnomcnological factor x~°- 57 of Eq. ([32| which makes the structure functions grow when x 
decreases and mimics in fact the total resummation 2 . In other words, the role of the higher- loop radiative corrections 
on the small- x behavior of the non-singlets is, actually, incorporated into SA phenomenologically, through the initial 
parton densities fits. It means that the singular factors can be dropped from such fits when the coefficient functions 
account for the total resummation of the leading logarithms and therefore fits for Sq become regular in x in this 

case. They also can be simplified. Indeed, if x in the regular part N (1 — x)@(l + "fx 5 ) of the fit |32|) is not large, 

all x -dependent terms can be neglected. So, instead of the rather complicated expression of Eq. (|32|) . 8q can be 
approximated by a constant or a linear form 



Sq(x) = iV(l + ax) 
with 2 phenomcnological parameters instead of 5 in Eq. (|32|) . 



(52) 



V. CORRECTING MISCONCEPTIONS 

The total resummation of ln fc (l/ir) allows to correct several misconceptions popular in the literature. We list and 
correct them below. 

Misconception 1: Impact of non- leading perturbative and non-perturbative contributions on the intercepts of g\ 
is large. 

Actually: Confronting our results and the estimates of the intercepts in Refs. [l8| obtained from fitting 
available experimental data manifests that the total contribution of non-leading perturbaive and non-perturbative 
contributions to the intercepts is very small, so the main impact on the intercepts is brought by the leading logarithms. 

Misconception 2: Intercepts of g\ should depend on Q 2 through the parametrization of the QCD coupling 
a s = a ( Q 2 ). 

Actually: This is groundless from the theoretical point of view and appears only if the the parametrization of the 
QCD coupling a s = a^k 2 ^) is kept in all ladder rungs. It is shown in Ref. [l(| that this parametrization cannot be 
used at small x and should be replaced by the parametrization of Eq. (I40|) . 

Misconception 3: Initial densities Sq(x) and Sg(x) are singular but they are defined at x not too small. Later, 
being convoluted with the coefficient functions, they become less singular. 

Actually: It is absolutely wrong: Eq. ([50")) proves that the pole singularity x~ a in the fits does not become weaker 
with the x-evolution. 

Misconception 4: Fits for the initial parton densities are complicated because they mimic unknown non- 
perturbative contributions. 

Actually: Our results demonstrate that the singular factors in the fits mimic the total resummation of ln fc (l/x) 
and can be dropped when the resummation is accounted for. In the regular part of the fits the x -dependence is 
essential for large x only, so impact of non-perturbative contributions is weak at the small- a; region. 

Misconception 5: Total resummations of ln fe (l/a;) may become of some importance at extremely small x but not 
for x available presently and in a forthcoming future. 

Actually: The efficiency of SA in the available small- a; range is based on exploiting the singular factors in the 
standard fits to mimic the resummations. So, the resummations have always been used in SA at small x in an 
inexplicit way, through the fits, but without being aware of it. 



2 We remind that our estimates for the intercepts AjvSi were confirmed (see Refs. [3) by analysis of the experimental data 
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VI. COMBINING THE TOTAL RESUMMATION AND DGLAP 



The total resummaton of leading logarithms of x considered in Sect. IV is essential at small-x. When x ~ 1, all 
terms ~ ln fc (l/:r) in the coefficient functions and anomalous dimensions cannot have a big impact compared to other 
terms. DGLAP accounts for those terms. It makes DGLAP be more precise at large x than our approach. So, 
there appears an obvious appeal to combine the DGLAP coefficient functions and anomalous dimensions with our 
expressions in order to obtain an approach equally good in the whole range of x : < x < 1. The prescription 
for such combining was suggested in Ref. [19| . Let us, for the sake of simplicity, consider here combining the total 
resummation and LO DGLAP. The generalization to NLO DGLAP can be done quite similarly. The prescription 
consists of the following points: 

Step A: Take Eqs. (f3"Tj) and replace a s by A of Eq. (|40|) . converting ^ns into jns and Cjfs mto Cns- 

Step B: Sum up the obtained expressions and Eqs. (|42|43|) : 

cns = C N ° S + Hs, Hns = Ins + H NS . (53) 

New expressions cns, h^s combine the total resummation and DGLAP but they obviously contain the double count- 
ing: some of the first-loop contributions are present both in Eqs. (f3"Tj) and in Eqs. (|42|43|l . To avoid the double 
counting, let us expend Eqs. (|42I43|) into series and retain in the series only the first loop contributions 3 : 



H 



(i) 

NS 



A{uoC h 
2V 



-1 

.U3 



C 



(i) 

NS 



A{luC f ) r 1 
2^ Vu? 



1 

2uj. 



(54) 



Finally, there is Step C: Subtract the first-loop expressions ([54]) from Eq. ([53]) ) to get the combined, or "synthetic" 
as we called them in Ref. fl9l . coefficient function cns and anomalous dimension hjys'- 



CNS = CNS 



C 



(1) 
NS' 



h NS = h 



NS 



H 



(i) 

NS- 



(55) 



Substituting Eqs. ([55]) in Eq. ([41]) leads to the expression for g^ s equally good at large and small x. This description 
does not require singular factors in the fits for the initial parton densities. An alternative approach for combining 
DLA expression for g\ was suggested in Ref. [20(. However, the parametrization of a s in this approach was simply 
borrowed from DGLAP, which makes this approach be unreliable at small x. 



VII. CONCLUSION 



We have briefly considered the essence of the IREE method together with examples of its application to different 
processes. They demonstrate that IREE is indeed the efficient and reliable instrument for all-orders calculations in 
QED, QCD and the Standard Model of EW interactions. As an example in favor of this point, let us just remind that 
there exist wrong expressions for the singlet oi in DLA obtained with an alternative technique and the exponentiation 
of EW double logarithms obtained in Ref. [1(| had previously been denied in several papers where other methods of 
all-order summations were used. 
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